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ABSTRACT 



Aims. We explore in this paper the ability to constrain the Galactic magnetic field intensity and spatial distribution 
with the incoming data from the Planck satellite experiment. 

Methods. We perform realistic simulations of the Planck observations at the polarized frequency bands from 30 to 
353 GHz for two all-sky surveys as expected for the nominal mission. These simulations include CMB, synchrotron 
\^ ■ and thermal dust Galactic emissions and instrumental noise. (Note that systematic effects are not considered in this 

paper). For the synchrotron and thermal dust Galactic emissions we use a coherent 3D model of the Galaxy describing 
its mater density and the magnetic field direction and intensity. We first simulate the synchrotron and dust emissions 
at 408 MHz and 545 GHz, respectively, and then we extrapolate them to the Planck frequency bands. 
Q ■ Results. We perform a likelihood analysis to compare the simulated data to a set of models obtained by varying the 

5— i ' pitch angle of the regular magnetic field spatial distribution, the relative amplitude of the turbulent magnetic field, 

^ | the radial scale of the electron and dust grain distributions, and the extrapolation spectral indices for the synchrotron 

, and thermal dust emissions. We are able to set tight constraints on all the parameters considered. We have also found 

that the observed spatial variations of the synchrotron and thermal dust spectral indices should not affect our ability 
to recover the other parameters of the model. 

Conclusions. From this, we conclude that the Planck satellite experiment can precisely measure the main properties 
of the Galactic magnetic field. An accurate reconstruction of the matter distribution would require on the one hand an 
improved modelling of the ISM and on the other hand to use extra data sets like rotation measurements of pulsars. 

"xl" ■ 
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Q| 1. Introduction 94 GH z. Furthermore. IBenoit et all (|2004l ) : iPonthieu et all 

fN| . (|2005l ) have observed significantly polarized thermal dust 

. The Planc k satellite (jPlanck- Collaboration! l2011at emission, up to 15 % in the 353 GHz Archeops channel. 

:• ; iTauber et all 120101. currently in flight, will provide The free- free emission is not intrinsically polarized and the 
. £h '• measurements of the CMB anisotropies both in tempera- anoma lous dust-correlated em ission is weakly polarized at 
ture and polarization over the full-sky at unprecedented 3+1-3 (jfiattistelli et all 12000) . At the Planck frequency 
$_i \ accuracy. It covers a large range of frequencies from bands the polarized contribution from compa ct and point 
C3 . 30 to 857 GHz and therefore is able to give a mea- sources is ex pected to be weak f or both radio (|Noltall2009ll 
surement of the foreground emissions. Planck makes and dusty ([Desert et all ll99Cf> sources. The frequency 
observations of the sky with a combined sensitiv- and spatial distributions of the polarized diffuse Galactic 
ity o f AT/Tcmb ~ 2 x 10~ 6 (jPlanck-Collaborationl emissions are not currently well known and the only avail- 
120051) and an angular resoluti on from 33 to 5 ar- able informations come from microwave and submillimeter 
cmin ([Planck- Collaboration! 120051 ) . In particular, because observations, 
of its 7 polarized channels it will for the first time allow 

the simultaneous precise measurement of the main sources The Galactic synchrotron emission originate from 
of polarized Galactic emissions: synchrotron and thermal re lativistic electrons spiraling along magnetic field lines, 
dust. The direction of polarization is orthogonal both to the 

line-of-sight and to the field lines. 

Usi ng the Wmap (W ilkinson Microwave Anisotropy 
Probe), IPage et all (l2007l l have shown that the synchrotron The synchrotron emission contributes principally to 
emission is highly polarized, up to 70% between 23 and diffuse emission at both radio and microwave observation 
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frequencies. Its spectral energy distribution (SED) is 
not known with accuracy however it is assumed to be 
well reproduced by a power law in antenna tempera- 
ture T v av^ B with a spectral index ranging f orm -2.7 to 
Hl987t iReich fe Reichl 119881: iReich et al ' 



3.4 (iLawson et al 



2004; iKogut et al 



20071 iGold et al.1 120091 



Fauvet et al 



20101 ) . The Galactic synchrotron emission has been well 



traced by the Leiden survey be t ween 408 MHz a nd 1. 4 
GHz (|Brouw fc Spoelstral fl976t iWolleben et all l200ll . 
the Parkes survey at 2.4 GHz dDuncan et al.1 Il999h and 
the MGLS ( Medium Ga lactic Latitude Survey) at 1.4 
GHz (jUvaniker et al.lll999D . At the microw ave frequencies 
it has been mapped by Wmap, see e.g. iHinshaw et al.1 
(f200^ ): |Pag"e et al.1 (l2007h . 

Thermal dust emission arises from dust grains in 
the interstellar medium (ISM) with typical sizes ~ 0.25 
/xm that are heated by stellar radiation (jPesert et al.1 
Il990lh This emission can be partially polarized as prolate 
dust grains align with their long axis perpe ndicular to 
the magnetic field (jDavis fe Greensteinl Il95l|) . The dust 
emission efficiency is greatest along the long axis, leading 
to partial linear polarization perpendicular to the magnetic 
field. The fractional polarization depends on the grain size 
distribution a nd is typically a few p e rcent at mill i meter 
wavelengths dHildebrand et al.1 119991: IVaillancourtl l2002t 
IFauvet et al.ll2010D . The thermal dust emission in intensity 
has already been well measu red by Iras form 5 to 100 
Aim (jNeueebauer et al.1 I1984D and COBE/FIRAS which 
provided the first polariz ed observation at high frequenc ies. 
Currently, Planck HFI (iPlanck-HFI-Core-TeamI [2010h is 
measu ring this emission in intensity ( Planck-Collaborationl 

i2oiibiicm . 

Based on the phys i cal ch aracteristics of the synchrotron 
emission, iPage et all (|2007l ) proposed a 3D model of the 
Galaxy including the distribution of relativistic electrons 
and the spatial distribution of the Ga lactic magnetic field. 
Independently, lHan et al.1 (|2004 120061 ) use d a 3D model of 
the f ree electrons density in the Galaxy (jCordes fe Laziol 
l2002t ) and a model of the Galactic magnetic field in- 
cluding regular and turbulent components to explain the 
observed rotation measu rements toward known pulsars. 
Based on these works ISun et aTl (|2008F ) performed a 
combined analysis of the rotation measurement of pulsars 
and of the p olarized W m ap d ata. This work has been 
extended by Jaffe et al I (I2010T) to study the Galactic 



plane and I Jansson et al.l (|2009T ) for the full sky. Recently 
IFauvet et al.1 (|2010[) proposed for the first time a coherent 
model of the synchrotron and thermal dust Galactic emis- 
sions using the Wmap and Archeops data. In addition 
to the above, many other related models and analyses can 
be found in the literature. 

We propose in this paper a method to study and con- 
strain the synchrotron and thermal dust polarized emis- 
sions using the Planck satellite observations that are cur- 
rently underway. Using realistic simulations of the Planck 
polarized data we forecast the expected constraints on a 
3D model of the Galactic magnetic field and the Galactic 
matter distribution. This paper is structured as follows: in 
Section l2~T1 we describe the models we used for the polarized 
components of the Galactic diffuse emissions and in Section 
3 we present the simulations of data used in this analysis. 



In Section 0] we describe our method to constrain these 
Galactic foreground emission models and we discuss the 
results in Section [5] Conclusions are presented in Section [BJ 



2. Models of polarized Galactic emissions 

A realistic model of synchrotron and thermal dust emis- 
sions can be constructed from a 3D description of the 
Galaxy including the ma tter distributi o n and the magnetic 
field structure. Following IFauvet et al.l (|2010fh we calculate 
the Stokes parameters I, Q and U of the emerging polarized 
Galactic emissions along the line-of-sight n. 
For the synchrotron emission, in galactocentric cylin- 
drical coordinates (r, 4 >, z) we use the following 
model (jFauvet et al.ll201C 



4ync(n) 

Qsync(n) 



0.408/ 



^Has/ff(n) 
^Has/ff(n) 

/ cos(2 7 (n, s))p s (Bf(n, s) + B t 2 (n, s)) n C R E (r, z)ds 



0.408 



E4ync(n) = / H as/ff(n) 



/ (Bf(n, z) + Bf(n, s)) n CRE (r, z)ds 
v, \A 



0.408 



/ sin(27(n, s))p s (B?(n, s) + ff 2 (n, sj) n C R. E (?-, z)ds 
f(B?(n, S ) + B?(n,s))n C RE(r,z)d S 

where B n (n,s) is the magnetic component along the line- 
of-sight n, and Bi(n, s) and B t (n, s) the magnetic field com- 
ponents on a plane perpendicular to the line-of-sight. The 
polar ization fraction p s is set to 75% ([Rvbicki fc Lightmanl 
1 1979th The polarization angle 7(11, s) is given by : 



7(n,s) 



arctan 



2B l (n,s)-B t {n,s) 
Bf(n,s) - £ 2 (n,.s) 



(1) 



The distribution of relativistic electrons, ncRE, is de- 
scribed in details in section 12.1.11 1 Has/ ff is a template 
temperature map obtained from the 408 MHz all-sky 
continuum survey (jHaslam et al.lll982lh This map is also 
included bremsstrahlung (free-free) emission. To substract 
this component we used the Wmap K-band free-free fore- 
ground map generated from the maximum entropy m ethod 
(MEM) (|Hinshaw et al.l 120071: iBennett et all l2003lh Note 
that t his template is not necessarily realistic: lAlves et al.l 
(|2010l ) have shown with radio recombination lines that in at 
least one region in the Galactic plane, this model appears 
to overestimate the amount of free-free. However, that will 
not have any impact in the following analysis because we 
used the same template in the simulated data and in the 
fitted models and the uncertainty at 408 MHz is very small 
compared to the synchrotron amplitude. The free-free map 
is then extrapolate d to 408 MHz a ssumin g a power law 
dependance as in iDickinson et~aTI (|2003l) . The spectral 
index /3 a used to extrapolate maps at various frequencies 
is a free parameter of the model and will be discussed later. 

For the thermal dust emission we used the following 
model (jFauvet et alJl2010h : 
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^dust(n) 
Qdust (n) 

^dust(n) 



I Ms (n) (|§) d , 
Ifds(n) (st)^ 

/ cos(27(n, s)) sin 2 (q)/ nol . m p rf n dust (r, z)ds 
J n dust (r,z)ds 

Ifds(n) (^)' d 

/ sin(27(n, s)) sin 2 (a)/ norm p d n dust (r, z)ds 



J n dust (r,z)ds 

whe re the dust p o lariza tion fraction p d is set to 10 
% (jPonthieu et al.l I2005D . and n dust (t, z) is the dust 
grain distribution discussed in section [2.1.11 The sin 2 (a) 
term accounts for geometrical supression and f n0 rm is 
an empirical factor which accounts for the misalign- 
men t between dust gra ins and the magnetic field lines 
(see iFauvet et al.l (120101) for d etails) . The reference map, 
Ifd s , is the iFinkbeiner et al. (I1999T) model 8 p rediction 
based on the Iras data ( Neugebauer et all 1 19841 ) and on 
the COBE/DIRBE data. The spectral index /3 d used to 
extrapolate maps at various frequencies is a free parameter 
of the model. This template seems to be a good repre- 
se ntation of the Archeops data at 353 GHz as discussed 
m M acias-Perez et al.l (|2007|) . Notice that for the polarized 
Planck frequencies, we are in the Rayleigh- Jeans domain 
and therefore, a power-law approximation for the dust 
intensity can be used. 



such that B = B reg + B tur b (bold letters indicate vec- 
torial quantities). We include only the isotropic part of 
the turbulent m agnetic field and no anisot ropic / ordered 
compo nent (see iJaffe et all (|2010h ). As in IFauvet et al.l 
(f2010h . our regular component is then equivalent to the 
sum of what Jaffe et al. (2010) call the coherent and 
ordered fields. For this analysis of synchrotron and dust 
emission only, the distinction is irrelevant. For the regular 



component we consider a Modified Lo 


2;arithmic Spiral 


model (MLS), presented in Fauvet et al.l ( 


201 0|) and based 


on the Wmap model (|Page et al.l 120071 
coordinates (r, <p,z) it reads : 


) . In cylindrical 



-Breg(r) 



£?(r)[cos(0 + 0) In I — sin(p) cos(x(r)) • u r 
Jo' 



- cos(0 + 0) In I — ) cos(p) cos(x(r)) • 
+sin(x(r)) ■ u z ], 



(4) 



where the pitch angle, p, is a free parameter of the model 
and /3 = 1/ tan(p). The scale radius tq is set to 7.1 kpc and 
x(r) = xo( r ){ z / z o) is the vertical scale height, with xo = 
22.4 degrees and z$ = lkpc. The intensity of the regular 
fie ld is fixed using p ulsar Faraday rotation measurements 
bv lHan et all (120061) : 



B{r) = B e 



2.1. 3D model of the Galaxy 

We describe here our 3D model of the Galaxy. 



where the large-scale field intensity at the Sun position is 
B Q = 2.1 ± O.ZfiG and R B = 8.5 ± 4.7kpc. The distance 
betw een the Sun and the G a lactic center, R Pl is set to 8 
kpc {EisenhaueretaLl (120031); iReid fc Brunthaier] (I2005D). 



2.1.1. Matter density 

We consider an exponential distribution of relativis- 

tic electrons npRF, on the Galactic disk motivated 

bv lDrimmel fc Spergell d200l : 



g °CRE,r 

ncRE^, z) = n , e -27— (2) 

COSh (Z/llCRE.h) 

where ncRE.r is the scale radius of the distribution and is 
a free parameter of the model. The vertical scale height, 
ncRE.h, is set to 1 kpc. The value o f no, e is set to 6.4 x 
10" 6 cm- 3 , following [SuneFaD (120081) . 

The distribution of dust grains ndust is described in the 
same way that of relativistic electrons with: 

e nd ' r 

ndust (r, z) = n , d r^T r, (3) 

cosh {z/n d ,h) 

where the scale radius nd_ r (r , z) is also a free parameter of 
the model (see IFauvet et al.l (|2010l ) for more details). The 
vertical scale height, n d , r , is set to 1 kpc. 

2.1.2. Galactic magnetic field model 

The Galactic magnetic field model is composed of two 
parts: a regular component and a turbulent component 



2.1.3. Turbulent component 

In addition to the large-scale Galactic magnetic field, 
Faraday rotation measurements on p ulsars in our vicinity 
have revealed a t urbulent component ([Lyne fc SmithllT989l : 
lHan et all |2004|) with an amplitude estimated to be of 
the same order of magnitude as that of the regular one 
(lHan et al.ll2006h . 

The turbulent magnetic field is assumed to be a 3D 
anisotropic gaussian random vectorial field and it is fully 
determined by a spherically symmetric power spectrum in 
the Fourier domain. Indeed, the magnetic energy Eg(k) 
associated with the turbulent com ponent can be describe d 
by a power spectrum of the form ([Han et al.l (|2004 12006|) ) 

E B (k) = c(£) a (6) 

where a = -0.37 and C = (6.8 ± 0.3) ■ 10" 13 erg cnT 3 kpc. 
More complex models of thi s anisotropic turbulent compo- 
nent have been proposed by (|Higdonlll98llCho et al.ll2002h 
and ICho fc Lazarianl (j2010l) but they are beyond the scope 
of this paper. Note also that we do not consider here the 
so-called ordered turbulen t component of th e Galactic mag- 
netic field as discussed in I Jaffe et al.l (j2010l ). 
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3. Simulated data 

We have performed simulations of the Planck polar- 
ized observations. We considered all the polarized chan- 
nels f or the LFI (iBersanelli et all I20I0I: iMandolesi et all 
20101: iMenella et all 1201 ll) an d HFI (lLamarre et all 1201 a 



Planck-HFI-Core-TeanJHblOl) instruments at 30, 44, 70 



100, 143, 217 and 353 GHz for two full-sky surveys (14 
months). We did not use the total intensity maps in this 
analysis. For this set of simulations we gene rated full-sky 
maps in the HEALPix pixelisation scheme (jGorski et all 
120051) at N S ide = 128. The input templates are simply de- 
graded to that resolution using the HEALPix tools. 

At a given observation frequency, ^, we consider that 
the polarized observations (X v ) are the sum of the syn- 
chrotron (X£ ync ) and the thermal dust emissions (X£ ust ), 
where X is the Q or U Stokes parameter. We also add a 
CMB contribution (Xq MB ) and noise (X^) so that can 
finally write: X" = X» ust + X£ ync + Xl T + XY 



L CMB 



spectral index map is shown in the left panel of Figure [TJ 

To simulate the spatial variations of the spectral index 
of the thermal dust emission we assumed a random gaus- 
sian distribution with a mean value of 0d,mean = 1-4 or 
Pd,mean = 2.0 and a standard deviation of 0.3. The right 
panel of Figure [TJ shows the resulting thermal dust spec- 
tral index map. Table [JJ summarizes the resulting types of 
simulations that we carried out. 

3.2. CMB emission 

We produced maps at N s id e = 128 of the expected 
CMB si gnal at all the Planck frequencies. We used 
CAMB (jLewis et alJl2000l) to compute the CMB tempera- 
ture and polarization angular power spectr a for the Wmap 
ACD M best fit model as estimated by (|Komatsu et all 
I2009T ) . We also took into account gravitational lens- 
ing effects and assumed a tensor-scalar ratio, r, of 
0.1 (jEfstathiou et al.ll2009l) . 



3.1. Polarized Galactic emission components 

The polarized Galactic emissions are simulated using 
the 3D model described in Section 12.11 The parame- 
ters of the 3D model o f the Galaxy have been taken 
from iFauvet et all (120100 bv compa rison with the avail- 
able Archeops dBenoit et all 12001) and Wmap 5-years 
data (jHinshaw et al. 20091) . The pitch angle is set to 
-30 degrees and the radial widths of the distributions 
of dust grains n<j,r and ultrarelativist electrons ncRE,r 
are set to 3 kpc. We computed simulated data using a 
model of Galactic magnetic field without (Simu I) or 
with a turbulent component (Simu II). In the second 
case we set the amplitude of the turbulent component to 
Aturb = 0.25 x B(r). The simulations are extrapolated 
to each of the Planck frequencies assuming constant 
spectral indices in frequency and using two different 
configurations, cither spatially constant (Simu Cst) or 
variable (Simu Var) accross the sky. We have produced at 
least 10 simulations for each of the different configurations. 

In the Simu Cst case we set the value of the synchrotron 
spectral index to S = —3.0. Concerning the thermal dust 
emission we test both upper and lower expected limits 
of the spectral index. We tested 1.4 as a lower limit, 
obtained by comparison of a grey body law model to the 
Archeops and Iris data in the galactic plane. The upper 
lim it tested is 2, foll o wing r esults obtained at high latitude 
by iBoulaneer et all (|l996l ) using the FiRAS and Dirbe 
data. 

For the Simu Var case, the maps of S and (3d are 
built in the following way. For the synchrotro n emission we 
used t he map generated using a MCMC fit bv lKoeut et all 
(|2007h at 23 GHz applying some corrections. We cut all 
pixels where S < — 3.5 and S > — 2 .8, fo ll owing the re- 
sults i n iPage et all ([2007D : ISun et all (|2008| ): IFauvet et all 
(|2010h . and assumed for those pixels S = — 3.0 + 0. 1 * n, 
where n is a gaussian random normal distribution. Note 
that S is the typical average value found in the litera- 
ture. It is also important to clarify that this spectral index 
map was not derived from the MEM analysis of the Wmap 
data. However, it can be used as a fair representation of the 
variability of the spectral index. The resulting synchrotron 



3.3. Noise 

Planck noise maps for each of the frequency channels 
were computed using the mean sensitivity per pixel given 
in Table [2] We assumed isotropic random gaussian noise 
accross the sky. This is not actually the c ase for the real 
Planck scanning strategy (|Dupad 120051 ). However this 
should not have any impact on the presented results as 
we are mostly signal dominated. 

4. Method 



Following IFauvet et all (|2010h . in order to compare the 
models of Galactic polarized emissions to the Planck data 
simulations we computed Galactic profiles in polarization 
using the set of latitude bands (in degrees) [0,30], [30,60], 
[60,90], [90, 120], [120, 180], [180,270], [270,330], [330,360]. 
Note that the intensity profiles are not used in this paper as 
the intensity maps were constructed using fixed templates. 

Galactic latitude profiles for the diffuse Galactic 
polarized emission models were computed for a grid 
of models obtained by varying the pitch angle, p, the 
turbulent component amplitude, A tur b, the radial scale for 
the distribution of electrons, ncRE,r and of dust grains 
nd t r, and the spectral indices S and 0d- The latter were 
assumed to be spatially constant accross the sky. All 
the others parameters were set to the values proposed in 
section 13.11 

We compared the simulated data sets to the Galactic 
emission models using a likelihood analysis where the 
log-likelihood function is given by: 



nfreq—l 



-logC 



Enjrt 
u=0 



E n^-i W („) - Ml x (n)f 



<x(nf 



•.turb 



(nf 



where the X are the Stokes parameters Q and U, and i 
and n index the longitude bands and the latitude bins, re- 
spectively. S" x (n) and M" x (n) are the set of simulations 
and models respectively for the polarization state X , i lon- 
gitude band and n longitude bin. The term a\ x (n) 2 is the 
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Table 1. Summary of the different types of simulations considered in this paper. 



Simulation 






Simu 1 


constant 


0.0 X B reg 


variable 


0.0 x B reg 


Simu 2 


constant 


0.25 X B reg 


variable 


0.25 x B r eg 




Fig. 1. Spatial variations of the spectral indices for the simulated synchrotron (left) and thermal dust Galactic polarized 
emissions (right). 

Table 2. Average la sensitivity per pixel (a square whose side is the FWHM extent of the beam) in ther modynamic 
temperature units, achievable after 2 full sky surveys (14 months) by Planck. (jPlanck-Collaborationl [20051 ) 



Center frequency [GHz] 


30 


44 


70 


100 


143 


217 


353 


(AT/T) polarization [fJ,K/K\ 


2.8 


3.9 


6.7 


4.0 


4.2 


9.8 


29.8 


Angular resolution [arcmin FWHM] 


33 


24 


14 


10 


7 


5 


5 



error associated to S" x (n) computed from the standard 
deviation of the data samples in each of the latitude bins. 
Note that it accounts both for the noise and signal disper- 
sion within the bin. The term cr^S"" (n) accounts for the 
additional variance due to the turbulent component of the 
magnetic field. Indeed, as the magnetic field is considered 
to be a random distribution, we need to take into account 
in the likelihood function an extra correlation matrix. We 
approximated this matrix to a diagonal one. We used 10 
simulations of the Galactic turbulent contribution at each 
Planck frequency band to estimate a^'^ rb (n). Note that 
the latter term is proportional to A tur b and also to the 
extrapolation term, (-[r~j)^ , both for the synchrotron and 
thermal dust components. This may introduce a small bias 
in the A tur b, Ps and (3d parameters. Increased turbulence 
can be balanced by a steeper spectral index and therefore 
it is possible to lower the resulting \ 2 f° r some parameter 
combinations. We will see, however, that the impact of this 
bias on the results remains small. 



5. Results and discussion 

The expected constraints on the parameters of the polarized 
Galactic emission models using the simulated Planck data 



are given in Tablcs@]and[5]for the various sets of simulations 
considered. 



5.1. Simu I 

The constraints obtained for the Simu I cases are pre- 
sented in the first line of Tables [4] and [5] The associated 
marginalized likelihood in 1 and 2D are shown in Figure [5] 
for the parameters A tur b, p, ncRE,r, n<i,r, (3 S and [3d, and 
we present the 1, 2 and 3c confidence level contours. We 
can see that there is no correlation between the parameters. 
We are able to tightly constrain all the parameters of the 
Galactic emission models. Furthermore, we can see that the 
best fit values coincide with the parameters used in the in- 
put simulations. Therefore, there is no indication of bias in 
the method. We can see that spatial variations of the spec- 
tral indices do not disturb the constraints on the parame- 
ters. Indeed the expected constraints on all the parameters 
of the models, including the values of spectral indices, are 
unchanged. This could be explained by the fact that we 
use a Galactic profile based comparison which is not very 
sensitive to pixel-to-pixel variations but rather to global 
features. The constraint on the dust grain density parame- 
ter is weaker than that for the relativistic electrons density. 
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Table 3. Parameters of the 3D Galactic polarized diffuse emission models. 



Parameters 


Range 


Binning 


P (deg) 


[-80.0,80.0] 


10.0 


-A-turb 


[0, 2.5] * B reg 


0.125 


ncRE,r (kpc) 


[1.0,20.0] 


1 


n d , r (kpc) 


[1.0,20.0] 


1 


Ps 


[-4.3, -2.4] 


0.1 


Pd 


[2.0,4.0] 


0.1 



Table 4. Best-fit parameters for the Galactic polarized emission models in the case of /3d = 1-4 in the simulated data. 
In the case of simulations with turbulence, Simu II, we give for each parameter the variance among the set of simulation 
results. 



Simulation 


Simu I 


Simu II 


ftsimu 


Cst 


Var 


Cst 


Var 


-^-turb 


< 0.1 


< 0.1 


0.25±^ 


0.25+S"? 


p(deg) 


-301* 


-3011 


-30l? 4 (5) 


-30t; 4 (5) 


rlCRE,r 


3+i- b 


3+i b 


121° (4.5) 


121° (4-5) 




3+} u 


3+{ u 


< 16 (5.1) 


12+6_ 8 (4.5) 


Ps 


n n +U.(Jb 

°- u -o.i 


o n +u.U5 

°- u -0.1 


-3.ll£ 7 ( < 0.05) 


-3.1+^ ( < 0.05) 


Pd 


1.4 ±0.1 


■■■•^-0.4 


1.5tS;§ (0.05) 


1.51^ (0.05) 



Table 5. Best-fit parameters for the Galactic polarized emission models in the case of /3d = 2.0 in the simulated data. 
In the case of simulations with turbulence, Simu II, we give for each parameter the variance among the set of simulation 
results. 



Simulation 


Simu I 


Simu II 


ftsimu 


Cst 


Var 


Cst 


Var 


-^turb 


< 0.1 


< 0.1 


0.35l'o i (0.06) 


0.351^(0.06) 


p(deg) 


-301* 


-30l| 


-301}* (7.5) 


-301^ (7.5) 


n,CRE,r 


3+{- b 


3+}- b 


121° (4.1) 


121° (4.1) 


n d , r 


3+{ u 


3+} l) 


< 16 (2) 


121° (2) 


Ps 


o n +u.U5 

°- u -0.1 


q n +0.U5 

°- u -o.i 


-3.11^7 ( 0.05) 


-3.11^7 (0.05) 


Pd 


2.0 ±0.1 


z - u -0.4 


2.2^:S (0.05) 


2.3t£S (0.05) 



This is probably due to the fact that the synchrotron is 
dominant at the Planck frequencies. 

5.2. Simu II 

The results concerning the Simu II cases, i.e. those includ- 
ing a turbulent component, are summarized in the lines 3 
and 4 of Tables 0] and [5l The marginalized likelihoods in 1 
and 2D for the parameters A tur b, p, ncRE,r, nd,r, B s and 
Pd are shown in figure |3] on which we present the 1, 2 and 
3(7 confidence level contours. We observe that there is no 
correlation between parameters, however, the constraints 
on the different parameters are weaker by a factor of 2 or 
more compared to the Simu I case. We observe a small 
bias on the best-fit values of the two spectral indices, but 
it is within the ler error bars; it is related to the additional 
noise-like term added to the likelihood calculation by the 
turbulent component of the magnetic field. Furthermore, 
the radial scales ncnE,r and rid, r are not constrained. A 
degeneracy between the matter distribution and the other 
parameters of the models is induced by the method of con- 
struction but is not visible in the Figure [3] Nevertheless, 
an upper limit can be set. Because the millimeter and sub- 



millimeter data are not very se nsitive to changes in those 
parameters, as was discussed in iFauvet et al.l (|2010D . only 
an upper limit can be determined. 



6. Summary and Conclusions 

We proposed a method to estimate the expected con- 
straints on the Galactic diffuse polarized emissions and 
the Galactic magnetic field at large scales using the 
Planck data. With this aim, we computed realistic 
simulations of the Planck data at the polarized fre- 
quency bands for two all-sky surveys. These simulations 
include CMB, synchrotron and thermal dust emissions 
and instrumental noise. For the synchrotron and thermal 
dust Galactic emissions we used a coherent 3D model of 
the Galaxy describing the magnetic field direction and 
intensity and the distribution of matter. The relativistic 
electron and dust grain densities were modeled using 
exponential distributions in galactocentric coordinates. 
For the Galactic magnetic field we co nsidered the Modifie d 
Logarithmic Spiral model discussed in lFauvet et al.l (|2010t ). 
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Fig. 2. Marginalised likelihood in 1 and 2D for the parameters A tur b and p, ncRE,r , nd.r, Ps and (3d for the (Simu I 
Cst) case and we present the 1 (68.8%), 2 (95.4%) and 3cr (98%) confidence level contours. The values of the parameters 
included in the simulated data are set to A tur b=0, p = —30°, ncRE,r = n d.r = 3kpc, (3 S = —3.0, (3d = 1.4 respectively. 



We performed a likelihood analysis to compare the 
simulated Planck data to a set of models obtained 
by varying the pitch angle of the regular magnetic field 
spatial distribution, the relative amplitude of the turbulent 
magnetic field, the radial scale of the electron and dust 
grain distributions as well as the extrapolation indices of 
the synchrotron and thermal dust emissions. We are able 
to set accurate constraints on most of the parameters 
considered. We have also found that the observed spatial 
variations of the synchrotron and thermal dust spectral 
indices should not affect our ability to recover the other 
parameters of the model. The presence of a turbulent 
component of the Galactic magnetic field decreases the 
discriminatory power of the method for all parameters 
but only in the case of the radial scales of the relativistic 
electron and dust grain distributions does it prevent an 
useful measurement. The small degree of bias the results 
in simulations including the turbulent component should 
not strongly affect the results with the real Planck data 
in polarization, since it remains small compared to the 



uncertainties. 



We can conclude that using the Planck data we should 
be able to constrain simultaneously the parameters of the 
models of synchrotron and thermal dust emissions. In par- 
ticular, we expect to constraint the direction of the Galactic 
magnetic field at large scales and the relative contributions 
of the regular and the istotropic turbulent component of 
the Galactic magnetic field without using external datasets. 
With respect to the current analysis, the constraints on the 
dust grain density parameters could be improved using also 
the total intensity data at the Planck HFI channels, from 
100 to 857 GHz. More generally, a more precise reconstruc- 
tion of the matter distribution in the Galaxy would require 
on the one hand an improved modelling of the ISM and on 
the other hand extra data s et s like rotat i on m easurements 
of pulsars lHan et al.l (|2004f ): ISun et al.l (|2008l) . These ro- 
tation measurement data along with total intensity should 
also help to constrain the ordered turbulent Galactic mag- 



7 



L. Fauvet et al.: Galactic polarized foreground for PLANCK 





0.0 0.5 1.0 1.5 2.0 -SO-60-40-200 20 40 60 5 10 15 20 




0.0 0.5 1.0 1.5 2.0 -30-60-40-200 204050 5 10 15 5 (0 (5 20 



.51 -2.5 -2.51 -2.5 

-3.5 ~H -3.5 1 

-,o| 



0.5 1.0 1.5 2.0 



0.0 0.5 1.0 1.5 2.0 



-80-60-40-20 20 40 50 




-30-60-40-20 20 40 50 




-4,0 -3.5 -3.0 -2.5 



4 




-4.0 -3.5 -3.0 -2.5 



0.5 1.0 1.5 2.0 



Fig. 3. Marginalised likelihood in 1 and 2D for the parameters A tur b and p, ncRE,r , nd,r, Ps and j3d for the (Simu II 
Var) case and we present the 1 (68.8%), 2 (95.4%) and 3cr (98%) confidence level contours. The values of the parameters 
included in the simulated data are set to A tur b=0, p = —30°, ncRE,r = n d.r = 3kpc, /3 S = —3.0, /3 c i = 1.4 respectively. 



netic field (see Jaffe et al, 2010), which has not been con- 
sidered in the current work. 
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